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Thermal and charge transport in diffusive normal metal(DN) / insulator / s-, d- and p-wave 
superconductor junctions are studied based on the Usadel equation with the Nazarov's generalized 
boundary condition. We derive a general expression of the thermal conductance in unconven- 
tional superconducting junctions. Thermal conductance, electric conductance of junctions and their 
Lorentz ratio are calculated as a function of resistance in DN, the Thouless energy, magnetic scat- 
tering rate in DN and transparency of the insulating barrier. We also discuss transport properties 
for various orientation angles between the normal to the interface and the crystal axis of supercon- 
ductors. It is demonstrated that the proximity effect does not influence the thermal conductance 
while the midgap Andreev resonant states suppress it. Dependencies of the electrical and thermal 
conductance on temperature are sensitive to pairing symmetries and orientation angles. The re- 
sults imply a possibility to distinguish one pairing symmetry from another based on the results of 
experimental observations. 

I. INTRODUCTION 

Thermal and electrical conductances are basic properties of a metal. The abilities to carry heat and charge currents 
are related to each other. At low temperatures thermal conductivity n of a metal is linear in T, i.e., k oc T, while 
the electrical conductivity a approaches a constant. As a result, the Lorentz ratio becomes a universal constant: 
L = n/crT = 7r 2 /3e 2 . This characteristic feature is called the Wiedemann- Franz (WF) lawi and has been observed 
in various electron systems2*24i£i& which may be described by the Fermi liquid theory. A violation of the WF law 
implies the breakdown of the Fermi liquid description of the electronic states in a metal. For instance, in the normal 
state of high-T c cuprates the violation of the WF law suggests the non-Fermi liquid description^, A validity of the 
Fermi liquid picture in cuprates is an open question even now. 

The heat and charge currents are also important characteristics of a superconducting state. Thermal conductivity of 
a bulk superconductor was first discussed theoretically by Bardeen et al£. Large amount of work was done on thermal 
and electric transport in contacts between normal and superconducting phases (N/S junctions). In the pioneering 
work by Andree^ a new type of quasiparticle scattering at the N/S interface was discovered, the so-called Andreev 
reflection (AR), which crucially influences quasiparticle transport across the interface at subgap energies. The AR 
causes the exponential decay of a thermal conductance across the N/S interface with decreasing temperature 9 -, while 
it facilitates the transfer of an electric charge 10 . As shown by Blonder, Tinkham and Klapwijk (BTK)afi, the AR 
leads to the doubling of zero bias conductance across transparent N/S interface at low T. More recently, Bardas and 
Averinii and Devyatov et alM calculated heat current by generalizing the BTK modeliS for the electric transport in 
N/S junctions. The applicability of these theories is, however, limited to junctions in the clean limit. 

In most of practical N/S junctions normal metals are in the diffusive regime. Therefore the effect of impurity 
scattering on the transport properties has received a lot of attention. In diffusive normal metal / superconductor 
(DN/S) junctions, the diffusive motion of quasiparticles at a mesoscopic length scale around the interface strongly 
modifies the transport across a junction interface because of interference effects^. For example, the appearance of 
the zero bias conductance peak (ZBCP) in DN/S junctions is a direct consequence of the interference effect of a 
quasiparticle in DNMii£iMiiiLi2ii2i 2 ii 2 Ii2 2 i 2 £i2i. To study such interference effects, the quasiclassical Green's function 
£ neor j i 25 J 26 4 27 nag b een widely used because of its convenience and broad applicability. Based on this formalism, 
theory of charge transport in DN/S junctions was formulated by Volkov, Zaitsev and Klapwijk (VZK) 28 . By ap- 
plying the VZK theory, a number of authors investigated theoretically charge transport in various proximity struc- 
tureo 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 i 36 i 37 i 38 i 39 . This work was based on the boundary conditions for the Keldysh-Nambu Green's 
function at the DN/S interface derived by Kupriyanov and Lukichev (KL)a£ from Zaitsev's boundary conditional in 
the isotropic limit. The KL boundary conditions were recently extended by Nazarov within the circuit theory*^ and 
applied by Tanaka et alM to the study of charge transport in N/S junctions with arbitrary interface transparency 
in order to investigate more complex structures. The Nazarov's boundary conditions coincide with the KL boundary 
conditions when transmission coefficients are sufficiently low, while the BTK theoryafli is reproduced in the ballistic 
regime. It was shown ins 3 , that a zero bias conductance peak (ZBCP) in low transparent N/S junctions transforms to 
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a zero bias conductance dip (ZBCD) with increasing the interface transparency. 

Heat transport within the quasiclassical approach was studied by several author si^£i2§i2L2£. In particular, Graf 
et a/4& calculated thermal conductivity for various unconventional superconductors. They showed that thermal 
conductivity of a clean two-dimensional d x 2_ y 2-wa,ve superconductor is proportional to T in the Born limit over a broad 
temperature range and proportional to T 3 in the unitary limit above some crossover temperature T* ~ 7, where 7 is the 
bandwidth of quasiparticle states bound to impurities. These results also suggest that thermal conductivity is sensitive 
to a pairing symmetry of an unconventional superconductor^. On the other hand, heat transport in N/S junctions 
has not received much attention so far, partially due to the lack of experimental data. Recently, sufficient progress 
was achieved which made it possible to study heat transport in unconventional superconductorsiSi^SiSiiSSiSi^iSi^ 
and this stimulated theoretical study of these phenomena. 

The formation of the midgap Andreev resonant state (MARS) drastically affects low energy transport in unconven- 
tional superconducting junctions££*SiSL§2i£i. It is now widely accepted that MARS is also responsible for ZBCP in 
N/S junctions of unconventional superconductors. To discuss interplay between the proximity effect and MARS, a new 
circuit theory for unconventional superconductors was proposed i n9 2 i9 3 i? 4 . In DN/d-wave superconductor junctions, 
MARS interfere destructively with the proximity effect in D N9 2 i9? . On the other hand, in DN/p-wave superconductor 
junctions the MARS and the proximity effect coexist*^. As a result, the conductance spectrum has a giant ZBCP 
and the local density of states in DN has a zero energy peak. We summarize relations between the proximity effect 
and the MARS in Fig. [TJ 
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FIG. 1: Schematic illustration of trajectories of incoming and outgoing quasiparticles at a DN/S interface with pair potential 
A±(<^). The pair potentials A± are given by A± = A(T) cos[2(0 =p a)] for d-wave superconductors (D-wave) where a denotes 
the angle between the normal to the interface and the crystal axis of d-wave superconductors. For p-wave superconductors 
(P-wave) the pair potentials A± are given by A± = ±A(T) cos[(0 =p a)], where a denotes the angle between the normal to the 
interface and the lobe direction of the p-wave pair potential. In the above, <j) denotes the injection angle of the quasiparticle 
measured from the :r-axis and A(T) is the maximum amplitude of the pair potential at a temperature T. 

The junctions in Fig. ^ are classified into four groups: (a) presence of the proximity effect, (b) presence of MARS, 
(c) absence of both MARS and the proximity effect, and (d) presence of both of them. Thermal and charge transport 
in N/S junctions can be described in terms of MARS and the proximity effect. 

The purpose of the present paper is to study the tunneling conductance, the thermal conductance and their Lorentz 
ratio in diffusive normal metal / insulator / s-, d- and p-wave superconductor junctions as a function of transparencies 
of insulating barriers at the interface, resistance Rd in DN, the magnetic scattering rate in DN, the Thouless energy 
Efh in DN, and orientation angles between the normal to the interface and the crystal axis of superconductors. 

The organization of this paper is as follows. In section II, we will provide the detailed derivation of the expression 
for the normalized thermal conductance. In section III, the results of calculations are presented for various types of 
junctions. They are applied to discriminate various pairing states. In section IV, the summary of the obtained results 
is given. In the present paper, we use the units with fcg = % = 1. 
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II. FORMULATION 



In this section, we explain the model and the formalism. We consider a junction consisting of normal and super- 
conducting reservoirs connected by a quasi-one-dimensional diffusive conductor with a length L much larger than 
the mean free path. The interface between the DN and the S has a resistance Rb while the DN/N interface has 
zero resistance. The positions of the DN/N interface and the DN/S interface are denoted as x = and x = L, 
respectively. According to the circuit theory, the interface between DN and S is subdivided into two isotropization 
zones in DN and S, two ballistic zones and a scattering zone. The sizes of the ballistic and scattering zones in the 
current direction are much shorter than the coherence length. The scattering zone is modeled as an infinitely narrow 
insulating barrier described by the delta function U(x) = H5(x — L). The transparency of the interface T n is given 
by T n = 4 cos 2 </)/(4cos 2 (j> + Z 2 ), where Z = 2H/vf is a dimensionless parameter, cf> is an injection angle measured 
from the interface normal and vf is the Fermi velocity. The interface resistance Rb is given by 

9 

Rb = Rq- 



I-w/2 d $ T ™ C0S 

where Rq is the Sharvin resistance R^ 1 — e 2 k 2 F S c j 'An 2 , kp is the Fermi wave length, and S c is the constriction area. 

We apply the quasiclassical Keldysh formalism in the following calculation of the tunneling and thermal conductance. 
The 4x4 Green's functions in DN and S are denoted by G\(x) and &2(x), respectively. The spatial dependence of 
Gi(x) in DN is determined by the static Usadel equation^, 

D ^P l {x)^^-] + i[H + it spin , G^x)] = (1) 

with the diffusion constant D in DN, where H is given by 

Ho 



H , 

Ho 

with H = ef3, Ti S pi n = ^T3Gi(x)t 3 is the self-energy for magnetic impurity scattering with the scattering rate 7 and 
e is the quasiparticle energy. The directions of magnetic moments of impurities are random. The self-energy is given 
by averaging with respect to directions of magnetic moments. Thus in our calculation G\(x) is unit matrix in spin 
space. 

The electric current is expressed using G\(x) as 

where (Gi{x)^§^-) K denotes the Keldysh component of {G x {x)^^-). 
The thermal current is also expressed as 



Hh 



L f°° deeTrKG^x)?^-)*}. (3) 



4e 2 i?d J dx 
It is convenient to use the standard 0-parametrization when function i?i (x) is expressed as 

Ri (x) — t 1 sin 0(x) cos ip + T2 sin 0{x) sin ip + T3 cos 6(x) (4) 

where cos^ = for singlet superconductors and sin?/; = for triplet superconductors (see Ref^). The parameter 
9{x) is a measure of the proximity effect in DN. 

Functions A\(x) and K\(x) are expressed as A\{x) = — t^R\{x)% and K\{x) = Ri(x)fi(x) — fi(x)Ai(x) with the 
distribution function f\{x) which is given by f\{x) = fi(x) + %ft{x). From the retarded or advanced component of 
the Usadel equation, the spatial dependence of 9(x) is determined by the following equation 

d 2 

D q^ ( x ) + 2 *( e + »7cos[0(s)]) sin[6(x)] = 0, (5) 



while from the Keldysh component we obtain 



D i [ ^ cos2Re[9{x)]] = - (6) 
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The average over injection angles of quasiparticles at the interface is defined as 

ctt/2 /.7r/2 



<B{cj>) >= f # cos <t>B(ct>)/ f ' 

J-tt/2 J -it 

with T{(f) = T n . 



d<j)T(cj)) cos cj) 

/2 J-rr/2 



A. s-wave case 

The Keldysh component K% (x) is given by K2{x) = R2(x)f2(x) — f2{x)A2(x) with the retarded component R2(x), 



the advanced component A2(x) and distribution function f2(x). Here R,2(x) = gf^ + fr 2 with g = e/y^e 2 — A 2 (T) 
and / = A(T)/y/A 2 (T) - e 2 , A 2 (x) = -t 3 rI(x)t 3 and f 2 {x) = f s = tanh[e/(2T)] in thermal equilibrium with 
temperature T. The boundary condition for Gi(x) at the DN/S interface is given by^, 

±(G 1 d -l±) lx = L _=R^<B>, (7) 
B _ 2T n [G 1 (L_),G 2 (L + )} 



4 + T n ([G 1 (L_),G 2 (L + )} + -2) 
For the electrical conductance, we obtain the following result at zero voltage^ 



1 f°° de 

M ) - £p / h 2r e U Eh I lid f L dx \ ^ > 

J0 C0Sn \-2T^<I ba > + L Jo cosh 2 8 irn (x)> 



<I b0 > L JO cosh 2 8 im (x) 

with 

r„ 2 A 1 + 2r»(2-r»)A 2 

h0 2 I (2-T n ) + T n \g cos 9 L +fsm6 L ] \ 2 ' 
Ai = (1+ I cos0 L | 2 + I sm6 L | 2 )(| g | 2 + I / | 2 +1) 



+4Imag[fg*]Imag[cos0 L sin0£], (9) 



A 2 = Real{g(cos L + cos 0£) + f(sin0 L + sin 0£)}, (10) 

where 9i m (x) and 0l denote the imaginary parts of 8(x) and 0(L ) respectively. 

Next we calculate the thermal conductance. Since DN is attached to the normal electrode at x = 0, 0(O)=O and 
fi(0) = fio with /;o = tanh[e/ (2(T + AT))] in thermal equilibrium at a temperature T + AT. 

The retarded part of Eq. Q is the same as in Ref. ]4Q The Keldysh part of Eq. J7J reads 

L ,dfi 2 < 7(,i > (fi(L-) - fs) . . 

^t(t^)cos Real[0(a;)] | x=i _= ^ , (11) 

with 

T 2 A' 1+ T„(2-T n )A 2 



4i = 



2 I (2 - T„) +T„[ 5 cos L + / sin 6 L ] | 2 ' 
A' 1 = (l+|cos0 L | 2 -|sin0 L | 2 )(| 5 |2-|/| 2 +1) 



+4Real[/ ff *]Real[cos 9 L sin 0£] , 



(12) 
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A 2 = 4Rea%] Real [cos 9 L ) - 4Imag[/]Imag[sin 9 L ] . (13) 
Substituting the results into Eq. @, we arrive at the final expression for the thermal current 

hk = L 2 r «(/■.-/"»* (») 



R b , R d r dx 



_k f 

(I bl ) L cos 2 Rc0(a;) 

and the thermal conductance 

K - fim r ^ fl 5) 

~ a™o AT " 2e 2 T 2 _/ / l V (15) 

r0 sh 2 f-M + ^ f ..^ 

cosn ^ 2T ; 1 ^ Jw j > l J cos -2 R c g( x ) I 
B. d-wave case 

In the following ± stands for the direction of motion along the x axis. We denote Keldysh-Nambu Green's function 
G2± as follows, 

where the Keldysh component Ki± is given by K%± = Ri±jiix) — f2{x)A2± with the retarded component i?2±, the 
advanced component Ai± and the distribution function f%{x). Here i?2± = ff±^3 + f±Vi with g-t = e/< /e 2 — Aj_(T), 



/± = A±(T)/y Aj_(T) — e 2 and ^2± = —^^±^3- The function /2(aj) is given in the previous subsection. For the 

electrical conductance, following RefS, we obtain the conductance, as(T), at zero voltage by replacing 1^2 into Ito 
in Eq. © where 

T T n Co 

J-b2 



2 I (2 - T n )(l + g + g_ + /+/_) + T„[cos0 L ( 5+ + 9 _) + sin0 L (/+ + /_)] | 2 

C = T„(l+ I cosfl L | 2 + I sm6 L | 2 )[| , 9+ + g_ | 2 + | /+ + /_ | 2 + | 1 + /+/_ + g + g_ | 2 + | /+«?_ - 5+ /_ | 2 ] 

+2(2 - T„)Real{(l + g* +g *_ + #/*)[(cos0 £ + cob 01) (0+ + <?-) + (sin0 L + sin^)(/ + + /_)]} 

+4TJmag(cos0 L sin^)Imag[(/ + + /_)(<£ + <?*)]. 
Let us calculate the thermal conductance. The boundary condition for G\(x) at the DN/S interface reads^ 



L . d6i 
R- d [Gl -dx- 



(Gi^) M _=i? ft 1 <I„>. (17) 



Here I n is given in the Appendix. The retarded part of Eq. I|17|) is the same as that in Ref. The Keldysh part of 
Eq. H17|) reads (see the Appendix) 

L ,dfi 2t} iivw mi <Ib3>{fl(L-)—fs) rio , 

R d ~(~dx~^ C ° S Rea ^°^ U=^-= > ( 18 ) 

with 

T r , C< 



r _ - » - 



2 | (2 - T„)(l + g+g_ + /+/_) + T„[cos0 L ( ff+ + <?_) + sin0 L (/ + + /_)] |* 
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For isotropic limit, where /+ = /_ and g + = g_ are satisfied, we obtain 1^ = hi. 

Applying a similar procedure as in the s-wave case, we obtain the thermal conductance by replacing h 3 into hi in 
Eq. |]2J 

For the ballistic limit, where 9l = Rd = 0, we can reproduce the formula for the thermal conductance for an ideal 
interface (T n = 1, i.e., Z = 0) found in the previous worki£: 

1 f°° e 2 < h3 > de 
K 2e*T*J i? b cosh 2 ^ 



where 



with 



A_ 



C. p-wave case 

Here, we restrict our attention to p-wave superconductors with S z — 0, where S z denotes the z component of the 
total spin of a Cooper pair. In this case the derivation is similar to that of the d-wave case. In the following, we will 
use the same notations as in the d-wave case. We can choose Ri — cos6(x)t3 + sin0(x)fi to satisfy the boundary 
condition at the interface. 2x2 matrix d {i = 1 — 6) is expressed in terms of linear combination of 1, fi, T2 and T3, 
and Bji = b\fi + &2T2 + ^373 with 

, _ ~T ln [T ln sm9 L + i(f+g- - f-g+)} 

°i 



(1 + I? n )(l + g+g _ + /+/_) + 2T ln [cos 9 L (g+ + gj) + 1 sm9 L (f+9- - /_<?+)] ' 

-T ln [T ln sinfl £ (l + g + . 9 _ + /+/_) + i(f+g_ - f_g+)] 

(1 + I* n )(l + g+g- + /+/_) + 2T ln [cos6 L (g + + gj) + zsin0 L (/+.g_ - /_<?+)] ' 

-T ln [T ln cosg £ (l + . g+ , 9 _ + /+/_) + g+ + g_] 

(1 + I? n )(l + g+g. + /+/_) + 2T ln [cos0 L ( ff+ + <?-) + i sm6 L (f + g. - /_<?+)] " 



(20) 



For the electrical conductance, following Ref. |6J, we obtain the conductance at zero voltage by replacing into 
ho in Eq. JBJl with 

j _ Cq_ 

b4 ~ 2 I (2-T n )(l+ 5+5 _+/ + /_)+T n [cos^( ff+ + 9 _)+isin^(/ +9 _- 5+ /_)] | 2 
C ' = T„(l+ I cos# L | 2 + I sin0 L | 2 )[| ,g+ + <?- | 2 + | /+ + /_ | 2 + | 1 + /+/- + g +9 - | 2 + | f+9- ~ 9+f- | 2 ] 
+2(2 - T„)Real{(l + g* + g*_ + f* + r_)[(cos6 L + cos6* L )(g + + 9-) + *(sin0 L + sin0 L )(/ + .g_ - <?+/_)]} 

+4T„Imag(cos 6> L sin6' i )Imag[z(/ + .g_ - g+f-){g* + + g*_ )]. 

Next we proceed with the discussion of the thermal conductance. The retarded part of Eq. i|17|) is the same as that 
in Ref. |^ We have to calculate the Keldysh part of Eq. I|17[) as in the d-wave case. Following equations are also 
satisfied: 

Dx (Ti„ - Ti n R~n + Tf n RiR^R p ) = B R , (21) 
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Br(1 + R m + TinRiRpR^) = Ti n R m Rp. (22) 

Applying similar procedure as in the d-wave case, we obtain the following expression for the thermal conductance 
by replacing Ib5 into Ibi in Eq. I|15[) where 

j _ ^0 

65 ~ 2 I (2 - T n ) (1 + .g+. 9 _ +/+/_)+ T n [cos L (g+ + g_ ) + i sin L (/+<?_ - |2 
C " - T n (l+ | cos^i | 2 - | sinfl L | 2 )(| g + + g_ | 2 - | /+ + /_ | 2 + | 1 + /+/_ + g + g_ | 2 - | - g + f_ | 2 ) 

+4(2-T n )[Real[(l + g* + g*_ + /+/*)(<?+ + 9 -)]Real(cos^)-Imag(si n ^)Imag[(l + g* + gl + r+f-Mf+9- - f-9+))]] 

+4T„Real(cos0 L sin0£)Real[i(/ +5 _ - f-g+){g*+ + g*)]. 

Taking into account all the above results, we can find simpler expressions for < Ibi > (t = 2 — 5) for any triplet 
superconductors (TS) and unconventional singlet superconductors (USS): 

< j > = < £o , . 

^ 2 | (2 - T n ) + T n (cos8 L g s + sm9 L fs) I 2 1 ' 



with 



for i — 2, 4 , or 



Do = T B (1+ | cos Ql I 2 + | sinfl L | 2 )[| 35 | 2 + | f s | 2 +1+ | f s | 2 ] 
+4(2 - T„)[Real(3s)Real(cos6iL) + Real(/ s )Real(sin0i,)] 

+4T n [Imag(cos0 i sm0£)Imag(/ s5s )] (24) 

D Q = T w (l+ I cos0 L | 2 - | sin0 £ | 2 )[| 35 | 2 - | f s | 2 +1- | f s | 2 ] 
+4(2 - T„) [Real(. 9s )Real(cos0 L ) - Imag(/ s )Imag(sin0 L )] 

+4T„[Real(cos 9 L sin0£)Real(/<jc/ s )] (25) 



for i — 3, 5 where 



„ _ > {9++. 9- )/(!+: 9+3- + /+/-) TS , , 

" 1 (9+ + ff-)/(l + .9+9- + /+/-) USS W 

r _ j *(/+.9- - /-.9+)/(l + 9+9- +/+/-) TS , , 

"'"^ (/+ + /-)/(! +.9+9- +/+/-) USS ^ 

f _ / (/+ + /-)/(! + 9+9- + /+/-) TS , , 

/S _ I Kf+9- 9+f-)/{l + 9+9- + /+/-) USS. ^ 

Negative sign appears in Eq. i|25[) . in contrast to Eq. I|24|l . since Cooper pairs cannot carry heat. 

In the following section, we will discuss the normalized conductance ctt(T) — ag{T) / ajq , the normalized thermal 
conductance kt(T) = k(T)/kn(T) and the normalized Lorentz ratio Lj- = kt(T) / '<jt{T) where cr/v and kn refer to 
the normal state and are given by ctat = l/(Rd + Rb) and kn(T) — tt 2 T/(3 (Rb + Rd) e 2 ) respectively. 
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III. RESULTS 



A. s-wave case 



First, we study the dependence of electrical conductance at the zero voltage, Ctt, on temperature as shown in 
Fig. |21 where we choose the relatively strong barrier Z = 3, Rd/ Rb = 1 hi (a) and (c), and Rd/Rb = 0.1 in (b) and 
(d). The Thouless energy is E Th /A(0) = 0.1 in (a) and (b), and E Th /A(0) = 0.01 in (c) and (d). The tunneling 
conductance has a peak at T/Tc — and minimum at T/Tc ~ 0.3 as shown in (a)-(d). The coherent AR due to the 
proximity effect in DN is responsible for the peak around zero temperature. The enhancement of the conductance 
is more pronounced for Rd/Rb = 1 in (a) and (c) than that for Rd/Rb = 0.1 in (b) and (d) because the proximity 
effect is more prominent for Rd/Rb = 1 than Rd/Rb = 0.1. The width of the peak is of the order of Erh/ A(0). The 
magnetic impurity scattering suppresses the proximity effect. Thus the height of peak decreases with the increase of 
7/A(0) as shown in Fig. [21 For low transparent interfaces the proximity effect enhances the conductance around zero 
temperature. For large T/Tc, <Jt increases monotonically with increasing T/Tc- 

The corresponding plots for Z = arc shown in Fig. [3] (a) - (d) where ctt has a dip at T/Tc — and a maximum 
at T/Tc ~ 0.5. It is known that similar dip-like structures appear also in the conductance as a function of a 
bias voltage^. The dip becomes broader for larger magnitudes of £r/i/A(0). For high transparent interfaces, the 
proximity effect suppresses the conductance around zero temperature. As a result, magnetic impurity scattering leads 
to enhancement of ctt, as illustrated by our numerical calculations. The non-monotonic temperature dependence of 
ctt is a unique feature of diffusive junctions. The structures in Fig.|31are essentially different from those by the VZK 
theory**, which stems from the high transparency at the interface. 
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FIG. 2: Normalized tunneling conductance for s-wave superconductor with Z — 3. (a) Rd/Rb = 1 and Eth/A(0) = 0.1. (b) 
Rd/Rb = 0.1 and E Th /A(0) = 0.1. (c) R d /R b = 1 and E Th /A(G) = 0.01. (d) R d /R b = 0.1 and E Th /A{Q) = 0.01. 

Next we study the thermal conductance, kt, as a function of temperature, where kt is normalized by its value 
in the normal state, and we fix Eth/A(0) — 0.1 and 7/A(0) = 0. We will show that kt is almost independent of 
£rfc/A(0) and 7/A(0), which implies that the coherent transmission due to the proximity effect does not affect the 
thermal conductance. Both a quasiparicle just on the Fermi energy and that with finite excitation energy (e) can 
carry the electric current. The results in Figs.[2]and 0show the sensitivity of electrical conductance around the zero 
temperature to Eth/ A(0) and 7/A(0) because the contribution of a quasiparticle just on the Fermi erengy is governed 
by the proximity effect. In the case of thermal conductance, on the other hand, only quasiparticles with finite energy 
can carry heat. At low temperatures, quasiparticles with e ~ T <C T c can contribute to kt- Such quasiparticles, 
however, are not allowed in the presence of the gap in superconductors. As a result, kt becomes almost zero around 
the zero temperature. In high temperatures such as T ~ Tc, only a quasiparticle with e ~ T c contributes to kt- In 
such energy range, the quasiparticle spectrum in DN is almost independent of Eth an d 7- Therefore kt is insensitive 
to Exh and 7. To substantiate it, we show some numerical examples in Fig. 0]where the independence is confirmed. 
As shown in Fig.El kt increases with increasing T/Tc for both Z = 3 and Z = 0. Contrary to the case of ctt, the 
magnitude of kt is reduced (enhanced) with the increase of Rd/Rb for Z — 3 (Z = 0). 
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FIG. 3: Normalized tunneling conductance for s-wave superconductor with Z = 0. (a) Rd/Rb = 1 and .Et;i/A(0) = 0.1. (b) 
Rd/Rb = 0.1 and E Th /A(0) = 0.1. (c) R d /R b = 1 and E Th /A{0) = 0.01. (d) R d /R b = 0.1 and E Th /A(G) = 0.01. 
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FIG. 4: (Color online) Normalized thermal conductance for s-wave superconductor with Z = 3 and Rd/Rb = 1- (a) 7/A(0) = 
and (b) E Th /A(0) = 0.1. 



We plot the Lorentz ratio, Lt, in Fig. |f)|for several Rd/Rb, where -&r/i/A(0) = 0.1 and 7/A(0) = 0. We confirmed 
that the Lorentz ratio is almost independent of i?Th/A(0) and 7/A(0). The Lorentz ratio is zero for small T/Tc an d 
is linear in T/Tc in the intermediate region. For Z = 3 and Rd/Rb = 0.1, Lt has a peak at T/Tc ~ 0-7, whereas it 
is a monotonic increasing function of T/Tc for Z = 3 and Rd/Rb = 1 (Fig- Et a ))- For Z = and Rd/Rb = L £t is 
linear in T/Tc for T/Tc > 0.3, while for Z = and Rd/Rb = 0.1 , is linear in T/T c in the intermediate region 

(Fig. EDO). 

Figure shows Rd/Rb dependence of tunneling conductance and thermal conductance normalized by their normal 
values at T/Tc = 0.8, where £t/j/A(0) = 0.1 and j/A(0) = 0. In low transparent interface (i.e., Z = 3), the 
probability of AR increases with increasing R d /Rb^. Thus the magnitude of ut increases with increasing Rd/Rb- 
On the other hand, kt is a decreasing function of Rd/Rb- In high transparent interface (i.e., Z = 0), the proximity 
effect suppresses the electrical conductance because the DN plays a role of the insulating barrier. In this case the 
probability of AR decreases with increasing Rd/Rb- Thus ctt is a decreasing function of Rd/Rb- The results also show 
that kt is an increasing function of Rd/Rb- In both (a) and (b), <jt and kt are close to constants independent of Z 
for sufficiently large Rd/Rb- 
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FIG. 5: Normalized thermal conductance for s-wave superconductor with i?Th/A(0) = 0.1 and 7/A(0) = 0. (a) Z — 3 and 
(b) Z = 0. 
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FIG. 6: Normalized Lorentz ratio for s-wave superconductor with _&rh/A(0) = 0.1 and 7/A(0) =0. (a) Z = 3 and (b) Z — 0. 



B. d-wave case 



In this subsection, we fix £r/,/A(0) = 0.1 and 7 = because kt and Lt are insensitive to these parameters. The 
pair potentials A±(T) are given by A±(T) = A(T) cos[2(0=F a)], where a denotes an angle between the normal to the 
interface and the crystal axis of d-wave superconductors, and cj> denotes an injection angle of a quasiparticle measured 
from the x-axis. The amplitude of pair potentials, A(T), is the same as that of s-wave superconductors. We choose 
< a < 7r/4. It is known that quasiparticles with 7r/4 — a < (f> < ir/4 + a can contribute to the MARS at the interface 
and are responsible for ZBCP in low transparent junctions. It was shown that the proximity effect and MARS do not 
coexist in the <i-wave symmetry. In fact, at a = 0, the MARS does not exist while the proximity effect is possible. 
On the other hand at a/ir = 0.25, the proximity effect is not possible, whereas the MARS appears (see Fig. ^ an d 
Ref. 63 ). Thus we can expect similar results to the s-wave symmetry in the case of a = 0. 

In Fig- El we plot the tunneling conductance as a function of temperatures for several choices of a, where Rd/Rb = 
0.1. At Z = 3 in (a), <tt for a/ir — 0.25 and 0.125 increase drastically with decreasing temperatures because of the 
resonant transmission through the MARS. Thus such behavior is not found in ar with a = 0. In high temperatures, 
<jt for a/ir = 0, 0.125 and 0.25 get close together. At Z — 0, ctt is a monotonic decreasing function of T/Tq for all 
a. The results show that ot slightly increases with a as shown in (b). The dependence on a of <jt at Z = is very 
small in comparison with that at Z = 3. 
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FIG. 7: Normalized tunneling conductance (a) and thermal conductance (b) for s-wave superconductor with T/Tc = 0.8, 
E Th /A(0) = 0.1 and 7/A(0) = 0. 
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FIG. 8: Normalized tunneling conductance for d-wave superconductor with Rd/Rb = 0.1, Eth/A(0) = 0.1 and 7/A(0) = 0. 
(a) Z = 3 and (b) Z = 0. 



In Fig. 13 we show thermal conductance as a function of temperatures for several a with Rd/Rb = 0.1. In both 
Z = 3 and 0, kt are monotonic increasing function of T/Tc and are proportional to T/Tc for small T/Tq- This 
linear dependence of kt for low temperatures is not seen in the s-wave symmetry (Fig. [SJ and reflects the line nodes 
of the pair potential. The results also show that kt depends on a for Z — 3 in (a), whereas it is almost independent 
of a for Z = in (b). 

In Fig. EH the Lorentz ratio is shown for several a, where Rd/Rb = 0.1. The Lorentz ratio has a peak at T/Tc ~ 0.4 
for Z = 3 and a = as shown in (a). This peak gradually disappears with increasing a. For Z — 0, Lt is proportional 
to T/Tc almost independent of a as shown in (b). 

Figures lTTl and lT^l show Rd/Rb dependence of tunneling conductance and thermal conductance which are normalized 
by their normal values at T/Tc — 0.8. For Z = 3 <tt is an increasing function of Rd/Rb for a/ir = as shown in Fig. 
1111 (a) because the probability of the AR increases with increasing Rd/Rb as in the s-wave symmetry. The peaks at 
Rd/Rb — in a/ir — 0.25 and 0.125 is a consequence of the MARS. The impurity scatterings in DN simply suppress 
the magnitude of <jt at a/ir — 0.25 because the the proximity effect is absent in this case. As a result, or becomes 
a decreasing function of Rd/Rb at a/n = 0.25 as shown in Fig. ^] (a). For sufficiently large Rd/Rb, &t is close 
to a constant. For a/ir = kt is a decreasing function of Rd/Rb while it is an increasing function of Rd/Rb for 
a/ir = 0.125 and a/ir = 0.25 as shown in Fig. ^2 (b). For sufficiently large Rd/Rb, kt is close to a constant for all 
a. In Fig. Illf bl. we find that the formation of the MARS suppresses Kt- This can be interpreted as follows. When 
the MARS is formed at the interface, the quasiparticle density of states becomes large around the zero energy. Such 
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FIG. 9: Normalized thermal conductance for d-wave superconductor with Rd/Rb = 0.1, -E Th /A(0) = 0.1 and y/A(0) = 0. (a) 
Z = 3 and (b) Z = 0. 
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FIG. 10: Normalized Lorenz ratio for d-wave superconductor with Rd/Rb = 0.1, J5yh/A(0) = 0.1 and 7/A(0) =0. (a) Z = 3 
and (b) Z = 0. 

states around the zero energy, however, can not carry the heat. The zero energy peak in the density of states means 
suppression of the density of states in higher energies which can carry the heat. Thus the formation of the MARS 
suppresses the thermal conductance. At Z = in Fig. 1121 the probability of the AR decreases as increasing Rd/Rb 
for all a. The line shapes of ay and Ky are understood in the same way as those in the s-wave case with Z = in 
Fig.H 

C. p-wave case 

Here we also fix _Ey/,/A(0) = 0.1 and 7 = because Ky and Ly are insensitive to these parameters. We choose 
the pair potentials A± as A± = ±A(T) cos[(c6 =F a)] where a denotes an angle between the normal to the interface 
and the lobe direction of the p-wave pair potential and A(T) is the maximum amplitude of the pair potential. In the 
following, we choose < a < ir/2. It is known that quasiparticles with injection angle cA with — 7r/2 + a < <f> < Tt/2 — a 
can contribute to the formation of the MARS at the interface. In particular at a = 0, the MARS and the proximity 
effect perfectly coexist, which causes the penetration of the resonant states into the DN§4. On the other hand, for 
a/ir = 0.5, neither the MARS nor the proximity effect exist (see Fig.^and Ref.|6J). 

In Fig. E| we show the calculated results of tunneling conductance for several a, where Rd/Rb = 0.1. At Z = 3, 
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FIG. 11: Normalized tunneling conductance (a) and thermal conductance (b) for d-wave superconductor with T/Tc = 0.8, 
Z = 3, E Th /A(0) = 0.1 and 7/A(0) = 0. 
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FIG. 12: Normalized tunneling conductance (a) and thermal conductance (b) for d-wave superconductor with T/Tc = 0.8, 
Z = 0, E Th /A(0) = 0.1 and t/A(0) =0. 

<jt for a/ir — and 0.25 increases drastically with the decrease of temperatures because of the resonant transmission 
via the MARS, while the result for a/ir — 0.5 monotonically increases with T. In the case of Z = 0, gt increases 
with decreasing T irrespective of a. 

The thermal conductance in Fig. ^]is a monotonic increasing function of T/Tc, where Rd/Rb = 0.1. Except for 
a = 0, kt is proportional to T for small T/Tc- This behavior is also found in the d-wave symmetry in Fig. [5] and 
stems from the line nodes in the pair potentials. At a = 0, kt is expected to be an exponential function of T/Tq. 
Thus kt increases with increasing a as shown in both Z — 3 and 0. Although line node of the pair potential exists 
in this case, kt has an exponential dependence on T as in the s-wave case because the direction of the line node is 
perpendicular to that of the thermal current. 

The Lorentz ratio has a peak at T/Tc ~ 0.3 for a — and Z = 3, as shown in Fig. I15f a~). This peak tends to 
disappear for larger a. At Z = Lt linearly increases with the increase of temperatures in intermediate temperature 
regime as shown in (b). 

Figures ITUI and IT7I display Rd/Rb dependence of tunneling conductance and thermal conductance at T/Tc = 0.8. 
The results are normalized by their normal values at T/Tc — 0.8. At Z = 3, <jt has a reentrant behavior in Rd/Rb 
for small a as shown in Fig. ^](a). It is noted that the finite energy states and the zero energy states contribute to 
the conductance in qualitatively different ways. In finite energies, ctt decreases with increasing Rd/Rb because the 
probability of the AR decreases with increasing Rd/Rb- On the other hand at the zero energy, <tt increase with the 
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FIG. 13: Normalized tunneling conductance for p-wave superconductor with Rd/Rb = 0.1, Eth/ A(0) = 0.1 and 7/A(0) = 0. 
(a) Z = 3 and (b) Z = 0. 
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FIG. 14: Normalized thermal conductance for p-wave superconductor with Rd/Rb = 0.1, Eth/A(Q) = 0.1 and j/A(0) = 0. 
(a) Z = 3 and (b) Z = 0. 

increase of Rd/Rb due to the formation of the resonant states^. For small (large) Rd/Rb, the contribution of the 
finite energies states (the zero energy states) dominates <jt- This explains the reentrant behavior in ax- In contrast 
to a — 0, <7t for a = 0.57T is almost constant as shown in Fig. ^| (a) because there are neither the proximity effect 
nor the MARS. The characteristics of the thermal conductance can be understood in a similar way. We note that 
the zero energy states never contribute to the thermal transport. At a = 0, kt increases with increasing of Rd/Rb 
in contrast to the d-wave case with a = 0. This difference stems from the existence of the MARS. At a/ir = 0.5, Kt 
is almost constant as shown in Fig. ^j](b). The line shapes of ctt and kt for a = and 0.25-7T at Z = in Fig. 1171 
are qualitatively similar to those for Z — 3. For a/ir = 0.5, or decreases with Rd/Rb and kt increases with Rd/Rb- 
These behaviors can be explained by the fact that the probability of the AR decreases as increasing Rd/Rb- From 
Fig. Ilfif b) for small Rd/Rb, we can also find that the MARS suppresses kt- 

On the basis of calculated results of electrical and thermal conductance, we propose a way to classify the pairing 
symmetries with several orientation angles into six groups as shown in Fig. 1181 We have studied seven junctions: 
s-wave superconductor (s), d-wave superconductor with a = (d(cv = 0)), d-wave superconductor with a/ir — 0.125 
(d(a/n = 0.125)), d-wave superconductor with a/n — 0.25 (d(a/7r = 0.25)), p-wave superconductor with a = 
(p(cv = 0)), p-wave superconductor with a/n = 0.25 (p(ev/7r = 0.25)) and p-wave superconductor with a/n = 0.5 
(p(cv/7r = 0.5)). We focus on the results with Z = 3 because symmetries of pair potentials are better characterized 
by transport properties in lower transparent junctions. Checking the line shapes of <jt and kt-, we can separate these 





15 



0.5 
1 

0.8 
0.6 

H 

1—1 0.4 

0.2 
0, 







/ a/ K-0 

o/tc-0.25 

; .' a/-=0.5 


(«) : 

Z = 3 


— — — f— > I i \ 

a/K-0 

al 71 =0.25 

o/tc-0.5 


' 1 1 


' (b) 






z=o 


0.2 0.4 0.6 


0.8 1 


T/T c 





FIG. 15: Normalized Lorenz ratio for p-wave superconductor with Rd/Rb = 0.1, -&rh/A(0) = 0.1 and 7/A(0) =0. (a) Z = 3 
and (b) Z = 0. 
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FIG. 16: Normalized tunneling conductance (a) and thermal conductance (b) for p-wave superconductor with T/Tc = 0.8, 
Z = 3, E Th /A(0) = 0.1 and 7/A(0) = 0. 



junctions into four groups as shown in the first and the second processes of Fig. I18f see Figs. |5J [5] |SJ El [E9 H4I and 
TABLE I). Next applying a weak magnetic field H parallel to the junction plane, ctt for d(a = 0) and p(a/7r = 0.25) 
junctions decreases as increasing magnetic field as shown in Fig. 1191 since the proximity effect is suppressed by the 
applied magnetic field 67 . On the other hand <tt for d(a/7r = 0.125), d(a/7r = 0.25) or p(ev/7r = 0.5) junctions 
is robust against applied magnetic field since there is almost no proximity effect (third process in Fig. I18fl . We 
note that the pair-breaking rate 7 is given by e 2 w 2 DH 2 /6, where w is the transverse size of the DN 6 ^. Assuming 
w = 10~ 5 m, D — 10 _3 m 2 /s, A(0) = 10~ 3 eV", and H = 10~ 2 T, we can estimate the pair-breaking rate 7/A(0) ~ 1. 
Distinguishing p(a/ir = 0.25) from d(a/7r = 0.125) is a delicate problem since both of them have the proximity effect 
and the MARS. The two junctions, however, have qualitative difference in the density of states (DOS) in DN. In 
p(a/7r = 0.25) junctions, the DOS has a zero energy peak due to the formation of the resonant states whereas the 
DOS in the d(a/7r = 0.125) junctions does not show such zero energy peak 6 ^. As a result, the conductance measurable 
by scanning tunneling spectroscopy (STS), ctq, reflects these features. Here ctq is defined as 



cr = 



1 

2T 



Re cos Ode 
cosh 2 [7^] 



(29) 



and Recos# is the DOS normalized by its normal states value. We plot it at x = 3L/4 as a function of T/Tc in Fig- 
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FIG. 17: Normalized tunneling conductance (a) and thermal conductance (b) for p-wave superconductor with T/Tc = 0.8, 
Z = 0, Er h /A(0) = 0.1 and 7/A(0) =0. 



TABLE I: Low temperature dependences of <7t and kt- 



pairing symmetry 






s 


reentrant 


exponential 


d(a = 0) 


linear 


linear 


d(a/7r = 0.125) 


inverse 


linear 


d(a/V = 0.25) 


inverse 


linear 


p(a = 0) 


inverse 


exponential 


p(q/tt = 0.25) 


inverse 


linear 


p(a/-K — 0.5) 


linear 


linear 



I2U1 In p(a/ir = 0.25) junctions, a peak appears at zero temperature in contrast to the case of the d(a/ir — 0.125) 
junctions. Thus we can easily distinguish these superconductors by STS. 
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FIG. 18: Chart for distinguishing s-, d- and p-wave superconductors. 



17 



d-wave 




FIG. 19: Magnetic field dependence of normalized tunneling conductance with Z = 3 and -&tti/A(0) =0.1 for (a) d-wave 
superconductor with a/ir = and (b) p-wave superconductor with a/ir = 0.25. 
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FIG. 20: The conductance measurable by STS at x = 3L/4 with Z = 3, i? d /i? 6 = 0.1, £ T h/A(0) = 0.1 and t/A(0) = for 
d-wave superconductor with a/ir = 0.125 and for p-wave superconductor with a/% = 0.25. 
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IV. CONCLUSIONS 



In the present paper, we have derived a general expression of the thermal conductance in normal metal / super- 
conductor junctions based on the Usadel equation under the generalized boundary condition. We have studied the 
electrical and thermal transport in diffusive normal metal / s-, d- and p- wave superconductor junctions in the presence 
of magnetic impurities in normal metals. The main conclusions are summarized as follows. 

1. The proximity effect does not influence the thermal conductance. This statement is illustrated with some 
numerically calculated examples. 

2. The midgap Andreev resonant states in d- or p-wave superconductor junctions suppress the thermal conductance. 
The formation of MARS drastically gathers the density of states at the Fermi energy near the DN/S interface. Such 
quasiparticles, however, do not carry heat because excitation energies of them are almost zero. 

3. The thermal conductance of the junctions reflects the existence of the line nodes of the pair potential except for 
the case that the direction of the line node is perpendicular to that of the thermal current. 

Electric conductance, thermal conductance, and their Lorentz ratio calculated as a function of temperature depend 
strongly on a pairing symmetry of a superconductor. This fact indicates a possibility of distinguishing one pairing 
symmetry from another by careful comparison of the present calculations and experimental results. 

In this paper, we have focused on N/S junctions of unconventional superconductors. So far an extension of the 
circuit theory to long diffusive S/N/S junctions has been performed by Bezuglyi et alp* 8 in s-wave symmetry. In 
S/N/S junctions, the multiple AR produces subharmonic gap structures in I-V curves^ ZWWWW. In S/N/S 
junctions of unconventional superconductors, it is known that MARS leads to the anomalous current-phase relation 
and temperature dependence of the Josephson current 77 . Effects of unconventional superconductivity of such I-V 
curves are an important future issue. The research in this direction is now in progress and the results will be reported 
elsewhere. 

The authors appreciate useful and fruitful discussions with J. Inouc, Yu. Nazarov and H. Itoh. This work was 
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V. APPENDIX 



Here we calculate the matrix current for the calculation of the thermal conductance of d-wave junctions. We denote 
H + , B n , I as follows, 
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h + = 7 y , h. = z n T > (so) 
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V o 


B A ) 




I A I 



B n — 

with H± = (G*2+ ± G*2-)/2. In singlet superconductors, we can choose R\ = cos 8(x)ts + sm6(x)f2 to satisfy the 
boundary condition at the interface. 

The matrix current is given by / = Tr„[/ n ] with /„ = 2[Gi, B n ] and 

B n = {-T ln [G ll HZ 1 ] + HZ 1 H+-Tf n GiHZ 1 H+Gi)-\T ln {l-HZ 1 ) (31) 

where 

T 



Tin — 



2-r n + 2VT=t: 

Next, we focus on the Keldysh component. We define lb 
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After straightforward calculations, 1^ is given by 

h = \Y, Tr ^ ~ Ai ^k + fi(L-)(Ri - M){B A - Br)]. 

n 

It is necessary to obtain Bk which is given by 

B K = D^N K - D^DkD^Na, (33) 

with 

D = -T^G^HZ 1 ] + HZ l H+ - T^G^H+Gu D = ( D * 



where Nk and N A is the Kcldysh and advanced part of N given by 



N = T ln -T ln HZ 1 +T? n G 1 HZ 1 H+ N 



N R N k 
N A 



We can express Nk and £>k as linear combination of distribution functions fs, fi(L_), and ft(L-) as follows, 

N K = difs + C 2 fi(L.) + C 3 f t (L_), 

D K = C 4 f s + <%/,(£_) + C 6 f t (L_), 

by 2 x 2 matrix Q (i = 1 - 6). Taking account of the fact that D^ 1 C 1 , b R x C 2l D^C^DZ^Na and D] i 1 C 5 D^ 1 N A 
can be expressed by the linear combination of f 2 and T3, while D~j^Cs and D^CqD^ 1 N A are proportional to the 
linear combination of 1 and fi, we can express as follows, 

h = ^Tr[fs(Ri - A 1 ){D R 1 C 1 - b^C i b- A l N A )+fi{L-){R 1 - A 1 )(D B 1 C 2 - D^C^Na) + (Ri - A 1 )(B A - B R )] 

n 

with 

C\ = Ti„L4~ - ii~ + T ln (R 1 R m 1 (R p - Ap) - R 1 (A m 1 - R^)A P )] 
C , 2 =T 1 2 „(i? 1 -i 1 )i m 1 i P 

C 4 = Ti„[iii(A~ - i?~ ) - (y^ 1 - ii~ )Ai] + Rm R P - AZnAp + T 1 2 „[-i?ii? m 1 (_R p - yl p )Ai + J?i(A~ - R^ n )A p Ai] 

C 5 = T ln [-(i?! - i^i- 1 + ^(fli - ii) - T^RiR^RpiRi - ii) + (4 - ii)^!" 1 !!)]. (34) 
We use the following equations 

Dx 1 (Ti n — Ti n R n i + T^RiR^Rp) — B R , 

Br(1 + #~ + TinRxRpR^) = Ti n R 71 ^Rp (35) 

where Br is given as Br = b\T\ + b 2 r 2 + 63T3 with 

-iT ln (f + g_ - /_#+) 



(1 + T?„)(l + g+g- + /+/_) + 27i„[cos0 L ( 5+ + + sin# L (/+ + /_ 
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-T ln [T ln sin0 L (l + g + g_ + /+/_) + /+ + /_] 



(1 + T?„)(l + g+g- + /+/_) + 2T 1 „[cos0 L ( 5+ + + sin0 L (/+ + /_)] ' 

, = -Tin Fin cos e L {\ + g+g _ + ;+;_) + g+ + g_] 

3 (l + T 1 2 n )(l+ 5+9 _ + / + /_) + 2r 1 „[co S L ( 5+ + ff _)+sin0 i (/ + + /_)]• 1 j 



Finally we reach the expression of given by the following equation. 

J =V— Cq(/s - h{L-)) 



C£ = T n (l+ I cos0 L | 2 - I sin0 L | 2 )(| ff+ +5- I 2 - I /+ + /- | 2 + I 1 + /+/- + 9+9- ? - I f+9- - 9+f- | 2 ) 

+4(2 - T„)[Real[(l + g* + g*_ + /;/*)(<?+ + . 9 -)]Real(cos# L ) - imag(sin0 L )lmag[(l + g* + g*_ + #/*)(/+ + /_)]] 

+4T„Real(cos L sin6>£)Real[(/+ + f-){g* + + g*_)}. 
This is a general expression applicable to any singlet superconductors without broken time reversal symmetry state. 
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